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Preroughening of close-packed fcc(lll) surfaces, found in rare gas solids, is an interesting, 
but poorly characterized phase transition. We introduce a restricted solid-on-solid model, 
named FCSOS, which describes it. Using mostly Monte Carlo, we study both statics, in- 
cluding critical behavior and scattering properties, and dynamics, including surface diffusion 
and growth. In antiphase scattering, it is shown that preroughening will generally show up 
at most as a dip. Surface growth is predicted to be continuous at preroughening, where 
surface self-diffusion should also drop. The physical mechanism leading to preroughening 
on rare gas surfaces is analysed, and identified in the step-step elastic repulsion. 



Thermal disordering in the height profile of a crystal 
surface will generally occur in two steps, as first shown 
by den Nijs. Using modified solid-on-solid (SOS) mod- 
els, JJ0] he found that regular roughening (temperature 
Tr) is preceded at some T c < Tr by a critical, prerough- 
ening (PR) transition. The intermediate disordered flat 
(DOF) phase presents a very special form of surface disor- 
der with proliferation of up-down correlated steps which 
causes the first layer to be only half-occupied. Based 
on the concept that PR is driven by latent tendencies 
towards reconstruction, (l]|| attempts at detecting PR 
focused mostly on metal surfaces, with uncertain results. 
|^,^) Surprisingly, the first convincing evidence of PR in 
unreconstructed surfaces came instead from rare gases. 
HQ The (111) faces of Ar, Kr, and Xe, exhibit re-entrant 
layering in adsorption isotherms above a T c « 0.85T m 
and below Tr ps 0.95T m (T m being the melting temper- 
ature), with a half-filled top layer between T c and Tr. 
This is a clear indication of PR, M and a simple hexago- 
nal model can in fact be constructed which contains PR 
with good similarities with the experiment. H 

However, the driving force for PR on Ar(lll) is cer- 
tainly not reconstruction, and remains to be understood. 
Moreover, the fcc(lll) surface, with its three sublat- 
tices, has a richer contents of degeneracy and symmetry 
than simple hexagonal, which needs to be addressed. At 
this stage, in fact, PR of fcc(lll) surfaces is not well 
characterized at all. We need, in particular, to under- 
stand scattering-related quantities such as order param- 
eters and susceptibilities, as well as important dynamic 
properties like surface diffusion and growth. This need is 
made more urgent by the strong possibility that PR (and 
not roughening) might more generally constitute the true 
onset of surface melting, as will be suggested by our sur- 
face diffusion results. 



We have constructed a restricted SOS model of 
fcc(lll) which yields well-defined answers to these ques- 
tions. We consider three triangular sublattices I — 0,1,2 
with ABC ABC . . . stacking. Heights in the l-th. sublat- 
tice are defined as integers hi = 3n.; + I, Vrij. Nearest- 
neighbor height differences are constrained to be ±1, ±2 
("infinite" bond strength). A positive energy is asso- 
ciated with height differences departing from the value 
attained in the perfectly ordered fcc(lll) surface (the 
ground state of the model). The hamiltonian is: 

H = J^2S(\hi - hj\ - 3) +K^5{\hi -h s \- 4) 

<2> <3> 

+ L^5(\h i -h j \-4) + M^8(\h i -h j \-5) •■• (1) 

where 53 (n) * s shorthand for sum over all pairs of n- 
th neighbors (distance measured in-plane), and the cou- 
plings J, K, L, M, . . . are positive energy parameters. We 
call this the FCSOS model (Fig. 1). Anticipating the re- 
sults, the minimal requirement in order to get a stable 
DOF phase into the model is to have at least non-zero 
K and L. Longer range interactions, including M and 
beyond, do not bring further changes. In this sense the 
FCSOS with L ^ is generic. 

We work out the essential features of the FCSOS phase 
diagram by the strip transfer-matrix method. Height 
parities (— l) hi are chosen to represent the SOS configu- 
rations of a N x oo strip. JhJ The fairly long-range in- 
teraction L obliges to compute all possible states of three 
spin rows. Our maximum strip size is N = 12. We lo- 
cate phase boundaries through the vanishing of the free 
energy of relevant interfaces. |Q Strictly at PR, we find 
that all off-plane excitations are costless. The step free 
energy vanishes at PR and remains zero inside the DOF 
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phase (due to proliferation of steps), whereas the cost of 
two parallel steps is non-zero until roughening. Within 
the transfer matrix framework these quantities are eval- 
uated by imposing suitable boundary conditions to the 
strip. Details of the calculation are reported elsewhere. 



We list a number of results (all the couplings from 
M onward are put equal to zero). When the only non- 
zero coupling in the model is J, there is no PR, and a 
simple roughening is found at e° J ~ 1.6. The situation 
remains unchanged when K > but L = 0. In particular, 
roughening takes place at e^ J ~ 1.5, when K = +oo. 
Finally, a PR transition is found along the K = L = +oo 
line at e^ J ~ 1.5. When L — K and J — the boundary 
between the DOF and the rough phase is found at e /3K « 
2. The overall phase diagram is sketched in Fig. 1. It 
resembles that of Ref. |lj , and is in fact more general (for 
positive couplings) as far as the nature of the phases it 
includes. 

Monte Carlo (MC) simulations confirm the above find- 
ings. We use N x N cells of increasing size, ranging from 
N = 24 to N = 96, and a standard Metropolis algo- 
rithm. After equilibration, more than two million MC 
sweeps are produced and average quantities such as the 
interfacial width (Sh 2 ), specific heat, and parity order 
parameter P = (3/iV 2 )(^ i (— l) hi ) , are evaluated. In 
an ideal DOF configuration the average height is half- 
integer, 5h 2 = 11/12 due to half-occupancy in the top- 
most layer, P should be zero, and its susceptibility \p 
should diverge at PR. 

As Fig. 2 shows, on the K = L = +oo line, P in- 
deed vanishes near e@ J = 1.5 and \p diverges. In the 
DOF phase, the average height is half- integer (Fig. 3), 
and correspondingly the occupancy of the top layer is 
one-half. Moreover, both the specific heat and its deriva- 
tive (not shown) remain finite at T c , thus signalling a 
higher-order transition than Ising, or 3-state Potts. The 
ordered phase and the DOF phase have the same degen- 
eracy (here 3-fold) . As a consequence, PR is expected to 
be non-universal j[| unlike Ising or 3-state Potts, where 
the degeneracy is fully removed above T c . 

Since the free energy of all step excitations is zero at 
PR, one expects the height fluctuations to diverge at T c , 
(Sh 2 ) ~ A ^ K InJV, but not below or above (until Tr). The 
MC results in Fig. 2 support this. Because of the Gaus- 
sian behavior at T c , critical exponents along the PR line 
can be recovered in terms of one parameter only, namely 
K c . At K = L = +oo and e 13 ' 1 = 1.5 we estimate 
K c « 1.06, whence v = ^ » 0.74, v = w 1.92, 

a = 2 - 2v w -1.85, (3 = « 0.71, and 7 = 

2 - a - 2/3 « 2.42. At finite L and K, the PR line 
is nearly insensitive to L/K, while roughening (where 
(Sh 2 ) ~ T7^ 2 \nN) shifts to higher (3K as L/K drops. 
Hence, the PR temperature is only controlled by the J 
value, which is then about 0.25 for argon (T c « 70 K), 



in units of the Lennard-Jones e (120 K for Ar). 

PR can be revealed in scattering experiments. 
]l|,||J|Jl^] Antiphase elastic X-ray or atom scattering is 
given by 

l(o,q z = £\ oc^q.g + ^xsCQ), (2) 

where a z is the vertical layer separation, Q and q z are the 
surface parallel and perpendicular momentum transfer, 
respectively. Here S = (3/N 2 )( J £2 i ai(—l) hi ) contains a 
shadowing factor on, |j) which may be taken, for example, 
to be one for local surface maxima and zero otherwise, 
and xs(Q) 1S tne local susceptibility. We find (Fig. 2) 
that the Bragg amplitude S 2 has a dip and vanishes at 
PR. However, unlike P, it is non-zero on both sides of 
the transition point (Fig. 2). This experimental signature 
can be used to detect PR on surfaces, including metals. 

Surface dynamics near PR is so far unexplored. We 
consider first surface growth. Continuous growth in the 
DOF phase is unlikely since the free energy of two paral- 
lel steps is non-zero. This implies existence of a free en- 
ergy barrier for the formation of a stable growth nucleus, 
and a behavior m oc e~ c / A ^ for the growth mobility, 
where Afi is the overpressure (smaller than a thresh- 
old A^ c ) and C is proportional to the square of the 
cost of two parallel steps. At finite size, we can also 
learn indirectly about the surface mobility m(N) from 
the behavior of the average surface height hit) at equi- 
librium. The surface as a whole diffuses in the form 
((h(t) - h(0)) 2 ) ~ 2d(N)t, with Einstein-like proportion- 
ality between d(N) and N~ 2 m(N). Fig. 3 shows both 
the equilibrium evolution of h (left), and the growth be- 
havior at finite A[i (right). At equilibrium the surface 
behaves like a Stokes particle undergoing Brownian mo- 
tion. The height is quantized, and is integer below T c and 
half-integer above T c . At PR quantization disappears, 
and d(N) is clearly much larger. This behavior is fully 
confirmed when A/i > 0. In both flat and DOF phases, 
growth is activated and occurs between quantized lev- 
els. At T c the surface is rough, and m soars accordingly, 
indicating continuous growth. 

The next dynamical issue is single-particle surface dif- 
fusion near the PR transition. On account of the criti- 
cal slowing-down affecting dynamical processes at a con- 
tinuous transition, and of a finite coupling of a migret- 
ing partucle to the DOF order parameter, we should 
expect a drop of the diffusion coefficient at PR. Jr| 
Dynamic scaling hypotesis combined with the assump- 
tion of Gaussian spreading of density inhomogeneities 
out of criticality predicts that diffusion should vanish as 

■y(*-2) 

D ~| t 2 -i , where z is a dynamical critical expo- 
nent and t = T Z Tc < 1. In the case of PR, the rugged 
landscape is expected to concur substantially to hinder 
particle hopping. The tracer diffusion D = (A? ,2 )/4t 
is extracted from a particle-conserving (Kawasaki) MC 
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simulation and displayed in Fig. 4. The size-dependence 
is maximum near T c , indicating a drop in the thermody- 
namic limit, as expected. A similar drop is found when 
approaching roughening from the DOF side (not shown). 

We are now ready to make contact with real fcc(lll) 
surfaces, such as Ar(lll). Modelling Ar by a Lennard- 
Jones 12-6 potential (a — tnn/1-0933, e = 120K, cut- 
off at r = 3.2cr), we consider a variety of step geome- 
tries, and use simulated annealing to optimize the T = 
geometry and energy. For both (100) and (111) single 
steps we find an energy (units of e) ~ 0.58 per unit 
lenght. This J value is about twice as large as that 
required for a PR transition at 70K. This is not sur- 
prising at all, in view of the additional vibrational and 
free-volume reduction factors which must be applied to 
T = step energies before using them near T m . We 
further find that two adjacent steps cost 1.14 (at the 
shortest distance) if up-down, and 1.49, 1.35, 1.22, . . ., 
2 • 0.58 + a/(l - l Q ) 2 (a/r& N ~ l-0> *o « -0.5), if up-up 
at distances 1/2, 1, 2, . . . , I lattice spacings, respectively. 
The lack of up-down interaction in FCSOS is thus jus- 
tified. The relatively short-range up-up repulsion needs 
to be qualified. To do this we determine effective FC- 
SOS K and L values (we take for simplicity K = L) so 
as to generate a probability distribution for the relative 
distance of two parallel steps which is closest to that gen- 
erated by the true long-range repulsion, using a method 
similar to Bartclt et al. fll4| For instance, when two par- 
allel steps are, on average, 4, 4.5, 5 lattice spacings apart, 
the long-range interaction a/(l — Iq) 2 enhances K and L 
from bare values of w 0.62 and « 0.16, up to larger effec- 
tive values K = L = 2.0, 14.5, > 10 4 , respectively. With 
reference to Fig. 1, we conclude that long-range elastic 
up-up step repulsion is crucial in giving rise to PR. This 
has a very transparent physical interpretation. Without 
step-step interactions (K = L = 0), PR and roughening 
initially coincide. Parallel-step repulsion leaves PR un- 
affected, but not roughening which is pushed at higher 
temperature. 

Now, some experimental consequences. The re-entrant 
behavior seen in rare-gas absorption isotherms, with half- 
filled layers above PR, agrees very well with our growth 
results (Fig. 3). Continuous growth strictly at T c is 
found, a remarkable feature of this system. We more- 
over find that antiphase scattering, if feasible, should 
show the critical drop at T c predicted for S 2 , followed 
by a recovery in the DOF phase, as in Fig. 2. This finite 
value of S 2 is generic, so long as the shadowing factor 
oti^X. Tracer diffusion along the surface, if measurable 
on a time scale short enough with respect to evapora- 
tion events, should also show a dip at T c , providing evi- 
dence for entanglement with critical fluctuations at PR. 
Finally, since parallel-step elastic repulsion is universal, 
we can expect that PR should also be common. In fact, 
we surmise that the abrupt vacancy proliferation gener- 



be precisely related to PR. It is hoped that these results 
will stimulate newer experimental efforts. 
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Figure Captions 



Fig. 1 Phase diagram of the FCSOS model (inset) for 
L = K. Thermal behavior of Ar(lll) corresponds 
qualitatively to the line K 3> J. 

Fig. 2 Preroughening critical behavior of the FCSOS 
model obtained by Monte Carlo along the line 
K = L = +00. From top to bottom: parity order 
parameter P, susceptibility \Pi Bragg scattering 
amplitude S, susceptibility xs, and interface width 
(5h 2 ) = (1/^V 2 )(E 4 (^-^) 2 )- Data for iV = 18(«), 
24(A), 36(D), 48(o), 60(*), 72(x), 96(solid A) are 
shown. Dotted lines are guides to the eye. 

Fig. 3 Average surface height h during a MC run, for 
N = 36, at three different temperatures along the 
K = L = +00 line. Left: equilibrium (A/z = 0). 
Right: growth (A/x = 0.0333). From top to bot- 
tom, eJ )J = 1.3 (DOF phase), 1.5 (PR), and 1.7 
(ordered flat phase). Note the continuous growth at 
PR, and the half-integer quantization in the DOF 
phase. Labels A, B, C describe the top layer (half- 
layer in DOF) sublattice. 

Fig. 4 Single-particle surface diffusion in FCSOS (K = 
L = +00). Data for N = 18(«), 36(D), 48(o), 
60(*), 96(solid A) are shown. Note the strong de- 
crease with increasing size, suggesting a dip at PR. 
The dotted line is a guide to the eye. 
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